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ABSTRACT 

We present a comprehensive nucleosynthesis study of the neutrino-driven wind in the aftermath of 
a binary neutron star merger. Our focus is the initial remnant phase when a massive central neutron 
star is present. Using tracers from a recent hydro dynamical simulation, we determine total masses 
and integrated abundances to characterize the composition of unbound matter. We find that the 
nucleosynthetic yields depend sensitively on both the life time of the massive neutron star and the 
polar angle. Matter in excess of up to 9 • becomes unbound until ^ 200 ms. Due to electron 

fractions of ~ 0.2 — 0.4 mainly nuclei with mass numbers A < 130 are synthesized, complementing 
the yields from the earlier dynamic ejecta. Mixing scenarios with these two types of ejecta can explain 
the abundance pattern in r-process enriched metal-poor stars. Additionally, we calculate heating rates 
for the decay of the freshly produced radioactive isotopes. The resulting light curve peaks in the blue 
band after about 4 h. Eurthermore, high opacities due to heavy r-process nuclei in the dynamic ejecta 
lead to a second peak in the infrared after 3 — 4 d. 

Subject headings: accretion, accretion disks — nuclear reactions, nucleosynthesis, abundances — neu¬ 
trinos — stars: neutron 


1. INTRODUGTION 

Since the first detection of neutron star binary sys¬ 
tems, the merger of two neutron stars or a neutron star 
with a black hole has been recognized as an ideal environ- 


sis ( 

Lattimer & Schramm 1974 

1 1976 Eichler et al.|l989) 

In the r-process, the neutron d 

ensity is so high that neu 


trons are captured much faster than the involved nuclei 
can decay. Once the neutron supply becomes insufficient, 
the temporarily formed, extremely neutron-rich nuclei 
decay toward the valley of beta-stability. Approximately 
half of all the elements heavier than iron are produced in 
this way. 

Matter that is ejected by hydrodynamical interaction 
and gravitational torques, often referred to as dynamic 
ejecta, has long been recogni zed as the likely source of 
the heaviest r-process nuclei (|Ereiburghaus et al.||1999b 

111 nnni H I 1/ ^1 


Rosswog et al.||1999| [Roberts et al.||2UII[ [Coriely et al 


vog 

2Q11[ IKorobkiii et al.||2Q12[ |B 


tokezaka et al.l 


et al.||2Qi31~[lIa 


|2013t|Rosswog|2013||Rosswog et al.|2014j^ 

box given nuclear physics input, thi s type of ejecta yields 


a very robust, “strong r-process” (S neden et al. 

dUanov et al. 


pattern due to efficient fission cycling _ _ 

Coriely et al.|2QlH[Korobkin et al.|2UI2[|Fetermann et al. 


2012 

1 

2015 
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et al.||2Q13[ |Wanajo et al.||2Q14| |Coriely 


Just et al. 2015a). This means that all merger 
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events produce practically identical abundance patterns 
(for this type of ejecta), independent of the details of 
the merging system. The abundance pattern, however, 
shows some sensitivity to the resulting distribution of fis- 
sion fragments and to the resulting beta-decay half lives 


Sion tragments and to tne resulting beta-decay nan liv e 
( Coriely et al.||2Q13 Eichler et al.||2Q15 Coriely 2015). 


With increasing sophistication of the included physics 
and of the numerical techniques also our view on the 
merger process and its mass loss channels has sharp¬ 
ened. Eor example, the dynamic ejecta have been rec¬ 
ognized to consist of a cold “tidal component” that is 
launched via gravitational torques and a hotter “inter¬ 
action component” that comes from the hydrodynamical 


interaction of the n eutron stars (Oechslin et al. 12007 


Ko- 


robkin et al. 

2012 

). Recent general relativity (CR) sim- 

ulations of tJ 

tie dynamic ejecta including neutrino trans- 

port ( 

Sekiguchi et al. 2015 Eoucart et al. 2015 

) show 

that shocks 

in the interaction component signihcantly 


increase the temperature. This hotter environment en¬ 
hances the weak reaction rates. In particular, the larger 
positron capture and electron neutrino absorption rates 
on neutrons increase the electron fraction, potentially up 
to values of 0.4. Therefore, the nucleosynthesis based on 
these general relativistic simulations can produce the full 
r-process from the fi rst peak (A = 80) to the third one 


{A = 195), see also Coriely et al. (2015). In contrast 


Newtonian simulations or approximate GR si mulations 


without weak reactions Coriely et al. (2011) 
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et al. (2013) only produce the heavy r-process elements 

In addition to the dynamic ejecta, the presence of two 
more mass loss channels has been appreciated and be¬ 
gun to be exp lored in more detail. These are neutrino 

driven winds (iRuffert et al. 1119971 iRosswog & Ramirez 
^ 1 I hW bAHUh ' ^ 1 I hWH 


RuizI 20021 IJJessart et al. |2o69[ T^rego et ^2014| |Jiist 
et al.l|2015ap and late-tinie accretion disk disintegration 


(iMetzger et ah||2QQ8l |BelobQrodov||2008| iLee et ai.||2009 

Fernandez lVletzger[2013f Just et al.j2015a ). The late¬ 


time disintegration of the accretion disk, driven by vis¬ 
cous heating and recombination of free nucleons into al¬ 
pha particles, is referred to as viscous ejecta in the fol¬ 
lowing. 

The historically favored production site, core-collapse 
supernovae, in contrast, seems to struggle to pro duce the 


conditions required for a successful r-process ( Arcones 


et al. 1120071 IHiidepohl et al.||2010t (Roberts et al.||2010t 

Fischer et al.M2010 

(Arcones fc Montes 201 1| Arcones 

fc Thielemannj2013 

). though some rare supernovae may 

plausibly contribute 

(iFujimoto et al.|2006, 20081 iWinteler 

et al.||2012 IVlosta et al.|2014 Nishimura et al.||2015|). 


sources of the heaviest r-process elements has recently 
received strong support from an observed near-infrared 
transient in the aftermath of a short gamma-ray burst 


et al 


(sGRB) (Tanvir et al. 2013 Berger et al 


.|2Q15[ ). This emission has be en interpreted as being 
a “kifonova” (Metzg er et al.|2Q ' lQ) (s ometimes referred to 
as “macronova”, cf. Kulkarni "(50051)), a transient that is 
powered by the radioactive decay of freshly synthesized r- 
process elements. In particular the luminosity peak after 
only several days and in the near-infrared - rather than 
after hours in the optical/UV as originally expected 0 
& Paczyhski 1998|) - supports the view that the ejecta are 
made of very heavy r-process elements that result in ef¬ 
fective opacities that are orders of magnitude larger than 
those of supernovae dKasen et al.||2013| Barnes fc Kasen 
2013 Tanaka fc Hotokezaka|20l3f|Grossman et al.|2014p 


IFis emission has be en interpre 


2013 Yang 


11 this interpretation is correct, this would lor the first 
time link compact binary mergers with GRBs observa- 
tionally, a long suspected, but actually so far unproven 
connection, and also with the formation of very heavy 
elements. Such kilonovae are an important facet in the 


multi-mes senger view of compact binary mergers ( Ross- ^ 
wog||2015 ) and, since they emit quasi-isotropically, they 


tional wave detections 

et al.|[2M3 Piran et al 


may become crucial in assuring the first direct gravita- 
d Metzg er fc Berger|2012||Nissanke 


wm 


Despite all of these promises, the role of compact bi¬ 
nary mergers for the chemical enrichment histories of 
galaxies, in particular at early times, is not yet suffi¬ 


cientl y well understood. While earlier work (Argast et al 
2004|) disfavored them as dominant source of r-process in 


the early galaxy, more r ecent studies are eithe r more opti¬ 
mistic about their role (Matteucci et al.|2014|) or actually 
favor them over super novae as major production sit e, see 
for example the recent hydro dynami c al stu dies of |Shen| 


et al. (2015) and van de Voort et aL] ( 2Q15| ). Here 


more 


work i s clea rl y indicated to settle th e case (see Gescutti 


Wehmeyer et al. ( 2Q15[ )). 


etal. (2015); 

Understanding the ditterent mass loss channels, their 
nucleosynthesis and possible electromagnetic transients 
is of prime importance for cosmic nucleosynthesis and 


galactical chemical evolution, for gravitational wave de¬ 
tection strategies and actually also for the GRB launch 
mechanism, since already a tiny mass loading of a fire¬ 


ball c an “choke” an emerging jet (Murguia-Berthier et al 
2014|). In this paper, we focus on the neutrino-driven 


wind from a post-merger remnant that consists of a cen¬ 
tral massive neutron star (MNS) surrounded by an ac¬ 
cretion disk and present a comprehensive nucleosynthesis 
study based on the fi rst three dimensiona l simulation re¬ 
cently performed by Perego et al. (2014). In the latter 
work preliminary nucleosynthesis results based on ten 
representative trajectories were presented. Here, we fol¬ 
low the nucleosynthetic evolution based on 17 000 trajec¬ 
tories and we discuss the implications for radioactively 
powered electromagnetic transients. 

The paper is structured in the following way. In Sect.[^ 
we briefly describe the hydro dynamical simulations, the 
tracers used to follow the ejecta and the physical input of 
the reaction network. We present the results of the nu¬ 
cleosynthesis calculations and report on time and angle 
dependency of the neutrino-driven wind in Sect. Fur¬ 
thermore, we consider mixing with the other channels of 
matter ejection and calculate the electromagnetic signal. 
Finally, we discuss our results and their implications in 
Sect. [4l 

2. METHODS 

2.1. Hydrodynamical simulations 

Our nucleosynthesis study is based on the first three di¬ 
mensional, Newtonian simulation of the neutrino-driven 
wind that emerge s during the afterma th of a binary neu¬ 
tron star merger (Perego et al.||2014). The merger rem¬ 
nant is characterized by a long-lived massive neutron 
star, surrounded by a quasi-Keplerian accretion disk with 
an inital mass of M^isk — 0.19Mq. For the present anal¬ 
ysis, we repeat and extend the previous simulation to 
longer times, and we follow a substantially larger num¬ 
ber of particles. 


We use the parallel grid code FISH (Kappeli et al. 
2011) to solve the hydrodynamical equations on a uni- 
torm Gartesian grid. Nuclear matter description is pro¬ 
vided by the TMl nuclear equation of state, supple- 
mented with electron-positron and photon cont ributions 
dTimmes fc Swesty||2000t |Hempel et al.||2012|). As ini- 
tial condition we use a late matter configuration from a 
three dimensional high resolution sim ulation of two non¬ 
spinni ng 1.4 Mq neutron stars (e.g. Price fc Rosswog 
(|2006|)). Neutrino-matter interactions are taken into ac- 
count using the multifiavor Advanced Spectral Leakage 
(ASL) scheme, see Perego et al. (2015, in preparation). 
It models the neutrino emission effectively by smoothly 
interpolating between diffusive and free-streaming rates, 
separately for different neutrino energies, and it has been 
carefully gauged at full transport calculations. Neutrino 
absorption is also included in optically thin conditions, 
based on the calculation of the neutrino densities out¬ 
side the neutrino last scattering surfaces. The list of the 
neutrino reactions impleme nted inside the ASL scheme 
can be found in table 1 of Perego et al. (2014). The 
present simulati on corresponds to an extension of the 
one presented in Perego et al. (2014). The simulations 
starts at 25 ms alter the beginning of the merger, in¬ 
cluding 10 ms of relaxation of the SPH final conditions 
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on the grid, and follows the post merger evolution for 
190 ms. The temporal evolution of the neutrino lumi¬ 
nosities and mean energies during the whole si mulation 
follows the same trends reported in Fig. 10 of |Pereg^ 


et al. (2014) for the first 90 ms. The relevant mean en¬ 
ergies stay approximately constant: ^11 MeV for the 
electron neutrinos and ^15.5 MeV for the electron an- 
tineutrinos. The total neutrino luminosities, integrated 
over the whole solid angle, decrease slowly with time. 
At the end of the simulation (tgim ~ 190 ms) the elec¬ 
tron (anti)neutrino luminosity due to cooling processes 
only is equal to ^ 2.2 • 10^^ erg/s 2.9 • 10^^ erg/s), 
while the inclusion of the neutrino absorption in opti¬ 
cally thin conditions decreases it to ^ 1.5 • 10^^ erg/s 
2.4 • 10^^ erg/s). Due to the larger optical depth 
along the equatorial plane, the neutrino fluxes measured 
far from the MNS along the polar direction are roughly 
3 times la rger than the ones i n the equatorial plane (cf. 
Fig. 12 of Perego et al. (2014|)). 

The inclusion of the energy deposition provided by the 
neutrino absorption on nucleons inside the disk drives 
a baryonic wind on a time scale of tens of milliseconds, 
see Fig. Matter expanding inside the wind becomes 
unboundat a distance < 600 km from the center. The 
resulting ejecta are confined within a polar angle of 60°, 
measured from the rotational axis of the disk. Due to the 
strong neutrino irradiation, the initially highly neutron 
rich matter inside the disk changes its electron fraction in 
the wind. The dominant z/e-a^bsorption of neutrons raises 


Yp tow ard equilibrium values (see e.g. Qian & Woosley 
( 1996| )). The evolution towards this equilibrium value 


may be affected by general relativistic treatment of the 
merger phase. GR simulati ons of a binary neu tron star 
merger including neutrinos (Neilsen et al. 2014 ) and fol¬ 
lowing the onset of the neutrino-driven wind ( [Sekiguchi 


et al. 2015) indicate that the luminosities can be larger 


than m Newtonian simulations, due to the higher tem¬ 
perature inside the massive neutron star and the disk. 
However, the ratio between the electron neutrino and an- 
tineutrino luminosities is similar in both cases, as well as 
the values of the mean energies. Therefore, the equilib¬ 
rium electron fraction is expected to be almost the same, 
while the evolution towards this Ye should be faster in 
the GR simulations. This difference does not come di¬ 
rectly from the Newtonian treatment of the disk and of 
the wind, but from the initial profiles resulting from the 
merger process. 

2.2. Tracers from the hydrodynamical simulation 

At the beginning of our simulation we place 10^ trac¬ 
ers inside our computational domain. The wind trac¬ 
ers are ejected from a region with initial density of 
P ^ lO^^gcm”^, as found i n backtracking pro cedures 
of preliminary tests following Perego et al. (2014). Thus, 
we locate our tracers initially inside the thick accretion 
disk where the density is between 2 • lO^gcrn”^ and 
2 • lO^^gcm”^ (see Fig. j^)- The number of tracers as¬ 
signed to each grid cell is proportional to its mass content 
and the actual location of each particle within a cell is 
randomly assigned. Hence, the resulting tracer distri¬ 
bution tracks the matter density distribution inside the 
disk. The mass of the region of the disk where the tracers 
are placed is 0.0542 Mq. Therefore, we assign to each 


of them an inherent mass of = m = 5.420 • IO^^Mq. 
Every tracer records the local properties of matter (den¬ 
sity, internal energy, electron fraction and velocity com¬ 
ponents) at its present location by linearly interpolating 
the corresponding values on the computational grid. The 
total specific energy, etot^ is computed as the algebraic 
sum of the kinetic, thermal and gravitational specific en¬ 
ergy. We recall here that the thermal energy is obtained 
from the relativistic internal energy, reduced by the rest 
mass energy (both obtained from the nuclear equation of 
state). All tracers are passively advected by the fluid and 
their location is evolved in time inside the grid by solv¬ 
ing the equation dx/dt = v with a second order accurate 
Euler integration scheme. A tracer particle is considered 
unbound if etot > 0 and if its radial component of the 
velocity is positive, steadily from a certain time until the 
end of the simulation. Eor each ejected tracer, neutrino 
fluxes and mean energies as a function of time come from 
the axisymmetric out put pro vided by the hydrody namic 
simulations, see Sect. 2.1 and Perego et al. (2014). 


2.3. Nucleosynthesis network 

Eor the post-processi ng of the tracers , a complete nu¬ 
cleosy nthesis network (Winteler 2012 Winteler et al. 
2012) is employed. Over 5800 nuclei between the val¬ 
ley of stability and the neutron drip line are considered, 
i.e. isotopes from H to Rg. The reaction rates are taken 
from the compilation of Rauscher & Thiele mann ( 200Q|) 
for the Einite Range Droplet IVlodel (4'RDlVl,|lVloller et al. 
( 1995|)). Weak interaction rates including neutrino ab¬ 
sorpt i on on nucleons are ta ken into account (Moller et al. 
20031 jErohlich et al. 2006). Eurthermore, neutron cap¬ 
ture for nuclei wi th Z > 80 and neu tron-induced fission 
rates are given by|Panov et al.|(|2010|) as well a s /d-delayed 
fission probabilities from Panov et al. (2005). 

We perform nucleosynthesis calculations lor over 17 000 
eject ed tr acers from the hydrodynamical simulation (see 
Sect. 2.2). Our computations start when the temperature 
drops below T = 10 GK. Then the initial composition 
is determined by nuclear statistical equilibrium (NSE) 
and is dominated by alpha particles, neutrons, and pro¬ 
tons. NSE is assumed to hold for T > 8 GK. Between 
10 GK > T > 8 GK, the network evolves the weak re¬ 
actions since they are not in equilibrium and change the 
Ye accordingly. As soon as the temperature undershoots 
the NSE threshold, the full network provides the abun¬ 
dances. The longest trajectories were simulated until 
4im ^ 200 ms, therefore w e extrapolate them foll owing 
the prescription outlined in Korobkin et al. (2012). 

Moreover, the energy generation by the r-process is 
calculated and its impact on the entropy is included 
(Ereiburghaus et al. 1999b). As is common practice, 
we neglect the feedback of nuclear heating on the den¬ 
sity evolution. In the case of the dynamic ejecta, it 
has been shown that this approximation is very good 
for nucleosynthesis calculations, although the long-term 
hydrod ynamical evolution o f the ejecta is heavily infiu- 


encec ^ [[Rosswog et al. 2014). The heating mainly orig¬ 


inates from /:/-decays and we assume that the energy is 

^ We expect that the impact on the nucleosynthesis in neutrino- 
driven winds is small, but since the expansion time scales here 
are different from the dynamic ejecta, this issue may need further 
scrutiny in the future. 
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Fig. 1.— X — z plane at 140 ms after the beginning of the simulation. The left panel shows the electron fraction with a superimposed 
contour plot of the density. A color gradient from white to black indicates the regions with densities of p = 10® g/cm^ to p = 10^^ g/cm®. 
Entropy profile and the projected velocity are presented in the right panel. The length of the arrows characterizes the magnitude of the 
velocity. 



Fig. 2. — Distribution of the tracers in the x — z plane at the 
beginning of the simulation. The density profile is shown for 
2 • 10® gcm“® < p < 2 • 10^® gcm“®. Tracers are marked by black 
dots, white dashed lines and numbers label the four angular re¬ 
gions of interest, which are referred to as (angular) bins 1 — 4 in 
the following. 

roughly equally distributed between thermalizing elec- 
trons and photons, a nd escaping neutrinos and photons 
( Metzger et aL|[2QlQ ). 

3. NUCLEOSYNTHESIS 
3.1. Time and angle dependeney 

The composition and amount of matter ejected in the 
neutrino-driven wind depends on the temporal evolution 
of the disk and the fate of the compact central object. 
Here we consider a long-lived MNS. In the following, we 
present abundances at various times after the merger 
and for different latitudes to understand the angular 
distribution and temporal evolution of ejected matter. 
The potential consequences for the mixing with the dy¬ 
namic ej ecta an d thu s for the light curve are discussed 
in Sects. [3^ and |3.3l respectively. In the following, tgim, 


refers to the time since merger. When abundances are 
shown at tgim, we consider all tracers having become un¬ 
bound until this time. Therefore, at later times early 
ejecta are also included. 

An idea of the amount of matter ejected and its poten¬ 
tial composition can be gained by exploring the depen¬ 
dence of the ejected mass on the electron fraction and 
time. This evolution is shown in Fig. Within the first 
50 ms only a marginal amount of mass with electron 
fraction < 0.3 gets unbound. Until tgim = 100 ms 
approximately 2 • 10 “^Mq are ejected with a central 
value of Yq ^ 0.3. After almost twice that time, i.e. 
4im = 190 ms, the mass of unbound material reaches 
9 • 10 “^Mq. These late-time ejecta has relatively high 
electron fractions (with a mean value between 0.3 and 
0.35) since neutrinos have more time to interact with nu¬ 
cleons. Furthermore, the growth rate of unbound mass 
by the late ejection of tracers remains approximately con¬ 
stant from 110 ms on. The substantial mass accretion 
rate inside the accretion disk, amounting to a few tenths 
of solar masses per second, provides quasi steady-state 
conditions with neutrinos as primary cooling agent. As 
a consequence, the neutrino luminosities decrea se only 
slowly over t he expansion time scale (cf. Fig. 10 in Perego 
et al. (2014)) and the interaction rates get close to equi- 
librium. 

In the following, abundances are shown using either of 
these two normalizations for a given nuclear mass A: 


Total mass: 


N 


X^ot{A) = J2MA) 


mi, 


( 1 ) 




with N, Xi and being the number of tracers 
considered, the mass fractions and the mass of one 
tracer i, respectively. This quantity has mass units 
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Fig. 3. — Evolution of the electron fraction for ejected tracers. 
The Ye distributions at various simulation times are represented 
by superposing bars with different colors. 


and thus allows to compare different ejection chan- 


nels to each other, as we will discuss in Sect. 3.2 


• Integrated abundances: 


Fint(A) = 


yM) ■ mi 

Ei=i rrii 


( 2 ) 


Normalizing the abundances via Eq. Q allows to 
explore the average nucleosynthesis yi^s of a cer¬ 
tain subclass of tracers, for instance tracers ejected 
into a particular solid angle. Additionally, the im¬ 
pact of single trajectories on the overall nucleosyn¬ 
thesis is directly evident from the comparison with 
the average abundance curve. 


Figure shows total ejected masses for simulation 
times tsim = 90 ms, 140 ms and 190 ms. The heavy 
nuclei beyond the second r-process peak {A ^ 130) are 
produced by early, very neutron-rich ejecta, as indicated 
by the overlap of the three curves at different simulation 
times. Later on, no heavy r-process elements {A > 130) 
are produced any more. There is a threshold Tg ^ 0.25 
below which heavy r-process elements (i.e., beyond the 
second peak at A = 130) can be synthesized (see also 


Kasen et al. (2015)). As time passes, a substantial 
amount of tracers contributing to nuclei with A < 120 


becomes unbound. This evolution is in fact a direct con¬ 
sequence of the trend of Yq presented in Fig. 

In addition, the nucleosynthesis has an angular depen¬ 
dency. We divide the neutrino-driven wind into four re¬ 
gions above and below the disk by setting cuts on the 
polar angle and investigate the angular dependence of 
the ejecta. Each of the cuts has a width of AO = 15°, 
hence the whole neutrino-driven wind is captured within 
0° < 0 < 60° (as a convention we define ^ = 0° at the 
poles). These four regions or bins are indicated in Fig. 
with white dashed lines and are labeled with the number 
of the bin. The properties of the tracers ejected until 
4im = 190 ms in each of the angular bins are shown in 
Figs. and when passing a sphere with a radius of 
750 km. We select the nucleosynthesis-relevant quanti¬ 
ties: electron fraction Tg, entropy per baryon s and radial 



Fig. 4.— Total final mass fractions times the mass of the ejecta 
until tgim = 90 ms, 140 ms and 190 ms. Early neutron-rich ejecta 
produce heavy r-nuclei with A > 130, leading to overlapping abun¬ 
dances in this region. Later ejecta contribute strongly to the lighter 
heavy elements (note the logarithmic ordinate). 

velocity v^. The major part of the cumulative mass is ap¬ 
proximately equally distributed to the two angular bins 
close to the accretion disk. On the contrary, the two bins 
in the polar region contain 15% — 20% of the total cumu¬ 
lative mass. While the entropy of the unbound tracers 
ranges from 10 /ce/baryon to 30 /cb/ baryon (Fig. [^, it is 
still very low and has therefore little impact on the nucle¬ 
osynthesis. The radial velocity provides a measure for the 
dynamical time-scale and is constrained to 0.04c — 0.08c 
(Fig.§. In the neutrino-driven wind the nucleosynthe¬ 
sis is most sensitive to the electron fraction distribution 
which varies for every angular region. As a general trend, 
the average electron fraction decreases as a function of 
the angle and reaches values down to 0.3 for the two 
zones closest to the disk. Both of these bins also con¬ 
tain extreme cases with very neutron-rich conditions of 

Fg < 0.2. 



Electron fraction Electron fraction 



Fig. 5. — Distribution of the unbound tracers in the Ye — s plane 
at a sphere with a radius of 750 km. We present the logarithmic 
scale of the total mass for all tracers in a certain angular bin of the 
Ye — s plane with a color code. The panels correspond to the angle 
intervals defined in Fig. ^ 

Figure shows the resulting integrated abundances for 
tracers ejected at tgim = 140 ms. All yields are presented 
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Fig. 6.— Same as Fig. but for the Ye — Vr plane. 

i n com parison to the solar abundances (dots, |Lodders 
(2003|)) for one angular bin in every panel. The thin gray 
lines represent individual tracers, while averages are in¬ 
dicated by a solid thick line in every panel. The patterns 
in Fig. [^reveal significant differences in abundances for 
distinct latitudes. We find that the first r-process peak 
{A ^ 80) forms for each angular region, whereas the sec¬ 
ond abundance peak (^4 ^ 130) is only reached in bin 3 
and bin 4. In particular, the angular zone closest to 
the disk, i.e., bin 4, successfully attains elemental abun¬ 
dances close to the solar system values up to the second 
abundance peak. When moving to lower polar angles 
(bin 1), less heavy elements are synthesized. 
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Fig. 7.— Individual abundances and integrated final abundances 
for ejecta until tgim = 140 ms in each of the four considered angular 
bins. Thick colored lines mark the average yields for each angular 
region. Note that we apply the same color code consistently in 
all following figures. Moreover, thin gray lines denote abundances 
from individual trajectories to hint the variety of the nucleosyn¬ 
thesis. Solar abundances are shown with dots for comparison. 


To gain further insights into the dependence on time, 
it is instructive to compare the trends of the angular 
bins at various times. Total masses as a function of mass 
number A (Eq. Q) for unbound tracers at 4im = 90 ms, 
140 ms and lOOnis are presented in Fig. At early 
times tsim ^ 90 ms, elements up to A ~ 120 are pro¬ 
duced predominantly by the ejecta in bin 3, while only 


bin 4 synthesizes nuclei within the vicinity of the third 
abundance peak {A ^ 195). The yields resulting from 
both bin 3 and bin 4 clearly dominate the abundances, 
as is expected from the distribution of the mass and the 
electron fraction for the tracers among the four bins (cf. 
Figs. and [^. In contrast, bin 1 and bin 2 give marginal 
contributions. With time, increasing yields for lighter 
heavy nuclei {A < 120) (cf. Fig. are present in all 
considered angular regions. 






A- 


bin 1: 0° 15° 

bin 3: 30° 45° 

bin 2: 15° 30° 

bin 4: 45° 60° 



A 


Fig. 8. — Total final masses. The panels show the yields for 
tracers ejected until tgim = 90 ms, 140 ms and 190 ms from top to 
bottom. We use the same color scheme as in Fig. □ to identify the 
different angular bins. 


Let us take a closer look at the relevant physical quan¬ 
tities characterizing the nucleosynthesis. Figure shows 
the evolution of individual electron fractions for tracers 
unbound at 140 ms (gray, thin lines), along with the cor¬ 
responding average curves (colored, thick lines). Note 
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that the elapsed time of a tracer is defined with respect 
to the moment when it gets out of NSE in the network 
calculations. The individual tracers provide information 
about the spread in the initial (i.e. at T = 10 GK), 
also visible in Figs. and The slight increase of 
the electron fraction Portly after the beginning of the 
calculations is due to neutrino absorption on nucleons, 
whereas the steep rise at around t^l0“^sis caused by 
/3-decays. 

The electron fraction in the neutrino-driven wind de¬ 
pends on the competition between the expansion and the 
weak equilibrium time scales, on the ratio between the 



the equilibrium Y^. Therefore, the very small amount 
of heavy elements (^4 > 130) that are produced in the 
initial phase of our Newtonian simulation may not be 
present in a GR simulation. In addition to the GR ef¬ 
fects, an accurate calculation of Y^ would require the 
usage of a detailed neutrino transport scheme, which is 
presently computationally prohibitive in long-time three 
dimensional simulations, without global symmetries. As¬ 
suming uncertainties of the order of 20% for the neutrino 
luminosities and of 10% for the mean energie^ we esti¬ 
mate a potential uncertainty of ^ 15% on the values of 
the equilibrium electron fraction. Given the broad range 
of Te obtained in the ejecta (Figs.|^[^and|^, we consider 
this uncertainty as being important to be mentioned, but 
not crucial for our analysis. 



Fig. 9. — Temporal evolution of the electron fraction for tracers 
that become unbound until tsim = 140 ms. The solid thick lines are 
averages (Ye) for every angular bin as introduced in Fig. Thin 
gray lines represent individual tracers. 


In order to examine the production of heavy elements, 
we present the evolution of neutron density the av¬ 
erage mass number {A) and the average proton number 
{Z) in Fig. 10 for different angular regions. At the be- 


^ This estimate is consistent with the quantitative differences 
spotted between our ASL scheme results and the ones obtained 
using a mult igroup flux limited d iffusion transport in axisymmetric 
simulations ( [Dessart et al.|2009| >. 


i30 3 


ginning, the neutron density is larger than 10 
and the composition is dominated by alpha particles and 
neutrons. The onset of the r-process nucleosynthesis is 
triggered when the temperature decreases to ^ 3 GK 
(marked by triangles in the top panel of Fig. 


M se 
nmctic 


e.g., the dependence of r-process efficiency as a function 
of entropy, electron fraction and expansion time scale. 


Hoffman et al.|[lM7| [Freiburghans et al.|[l999a| ). For 
t (1.5 — 3.0) • 10“^ s, most of the neutrons are con¬ 
sumed as indicated by the rapid neutron density drop. 
After this, matter beta-decays to stability and this leads 
to an increase of Yq (see Fig.[^. Then, the mass number 
stagnates, since no more heavy elements are formed (see 
middle panel of Fig. 10). No fission takes place, resulting 
in the monotonic evolution of mass and proton number. 
The higher neutron density in bin 4 favors the build-up 
of heavier nuclei compared to the other regions. With 
decreasing polar angle the neutrons are consumed ear¬ 
lier, as their initial density is smaller and the expansion 
velocity of the corresponding tracers is larger. 



Fig. 10. — Evolution of the neutron density Un (top), average 
mass number (A) (middle) and average proton number (Z) (bot¬ 
tom). The solid lines follow the mean values of the four angular 
bins (see also Fig. 1^. Triangles indicate times, when the temper¬ 
ature reaches T ^^GK, i.e. the onset of the r-process. 


3.2. Neutrino-driven wind and dynamie ejeeta 

The neutrino-driven wind is only one of three 
nucleosynthesis-relevant ejecta in neutron star mergers. 






































In order to have the complete nucleosynthesis picture, 
one would need to follow the evolution of the dynamic 
and disk (neutrino- and viscous-driven) ejecta, and in¬ 
vestigate how they mix. However, this is currently too 
complicated to be studied in a single simulation due to 
the different time-scales and physics involved. The dy¬ 
namic ejecta expand very fast from the beginning. While 
neutrinos may be initially important for the Yq their im- 
pact is insignificant when the disk becomes transparent 


( Fernandez fc Metzger |2Q 13 


& Fernandezl2014 


ust et al.J|2U15a[ ). JNevertheless, mixing of the diherent 


types of ejecta can occur once the unbound matter be¬ 
comes considerably decelerated by the ambient medium. 
Where this happens depends on the merger location with 
respect to the galaxy. Simple estimates indicate that 
the time scale for being noticeably decelerated is by far 
greater tha n years, see, for exam ple, the discussion in 
Sect. 2.1 of Rosswog et al. (2014|)). 

In the first stu dy of the mixing between the disk and 
dynamic ejecta, |Fernandez et al.| (|2015|) showed that 
there is no significant mixing between these ejecta. This 
study focused on neutron star black hole mergers, where 
neutrino-driven ejecta are less import a nt as it was dis- 
cussed also by Fernandez fc Metzger (2013) compared 


Metzger fc Fernandez ~ pOTT ). However, tneir conclu- 


to 

sion may also apply to our results that are based on a 
long-lived MNS. F rom the nucleosynthesis perspective. 
Just et al.| (|2015a|) have also investigated the mixing of 
disk and dynamic ejecta. They have modeled the long¬ 
time evolution of a black hole surrounded by a disk where 
matter is neutrino- and viscous-driven ejected. In a post¬ 
processing step, the nucleosynthesis of such simulations 
was added to the dynamic ejecta from a merger simula¬ 
tion ( Banswein et al.|[2Q13 ). 

We combine the nucl eosy nthesis of the neutrino-driven 
wind pres ented in Sect. |3.1| with t he dynamic ejecta com¬ 
puted by Korobkin et al.| (2012). Although the ejecta 
masses show a substantial spread depe nding on the pa - 
rameters of the merging binary system (|Rosswog||2013|), 
the resulting abunda nce patterns have been found to be 
practically identical ([Kor obkin et al.||2Q12). In contras t 
to previous studi es (i.e., kernandez & Metzger (2013); 
Just et al. (2015a)) with a black hole as central object, we 
assume a long-lived MNS. This has strong consequences 
on the amount and properties (i.e ., Te) of the matte r 
ejected by neutrinos as discu s sed by|Perego et al.j (|2014|) ; 
Metzger & Fernandez|(|2014|); |Kasen et al.|(|2015|. 

4'he comparison between dynamic and wind ejecta is 
given in Fig. 11 showing the total ejected masses (Eq. 
as a function of A for both contributions. Note that the 
approach here differs from the one in the previous sec¬ 
tions, in which we treated the different times tgim to be 
merely snapshots of the neutrino-driven wind. Now, we 
assume that the MNS collapses after 4im: terminating 
the ejection of further material in the wind. Three dif¬ 
ferent simulations times are considered, tgim = 90 ms, 
140 ms, and 190 ms. This comparison indicates that 
the wind ejecta complement the dynamic ejecta by pro¬ 
ducing elements below the second peak. We note that 
if the dynamic ejecta could produce the full r-process, 
as suggested by recent GR simulations (Wanajo et al. 


2014r[Sekiguchi et al. 
/I/ ^ _ A k I kvH -i- 


( iGoriely et al.|2015t , t 


2015) and by parametric studies 
le combination of the dynamic and 
wind nucleosynthesis would again allow the production of 


lighter heavy elements {A < 130) in addition to the heavy 
r-process. In the latter case, the wind ejecta would in¬ 
crease only the relative contribution of the lighter heavy 
elements. Moreover, the amount of wind ejecta becomes 
comparable to the dynamic ejecta if the MNS survives 
long enough (4im ^ 140 ms in our model). These ejecta 
amounts will be further enhanced by viscous ejecta that 
can pr oduce r-process ele ments from the first to the third 
peak (Just et al. 2015a). If the three ejecta completely 
mix, tne wind contribution may still lead to variations in 
the abundances below A = 130 (i.e., Z < 50). Such vari¬ 
ations are also expected from observations of the oldest 
stars (see Sneden et al. (2008) for a review). The size of 
the observed variations could help to constrain different 
contributions of the three ejecta types to the r-process 
abundances. 



Fig. 11.— Comparison of the nucleosynthesis yields for neutrino- 
driven wind (cf. Fig. and dynamic ejecta (labeled with dy¬ 
namic). While the dynamic ejecta produce very heavy nuclei, the 
neutrino-driven wind complements its abundances by producing 
the lighter heavy elements from the first to the second r-process 
peak. 


The oldest observed stars were formed in an interstel¬ 
lar medium polluted by only one or few nucleosynthesis 
events. Therefore, their atmospheres present a unique 
fingerprint for the r-process. Observations indicate that 
there are at least two types of abun dance patterns amon g 
the oldest r-process stars ( see, e.g., [Sneden et al.|(|2008|); 
Qian & Wasserburg (2007)): 1) stars with high enhance¬ 
ment of heavy elements [Z > 50) present a robust pat¬ 
tern for those and some variations for the lighter heavy 
elements {Z < 50); 2) stars with low enhancement of 
heavy r-process. Two stars are typically identified with 
these trends: 1) CS 22829-052 (Sneden et al.[[2003[) and 


2) HD 122563 ([Waller stein et al. 


2006). Recently 


Hansen et al. 


Tm 


i 


Honda et al. 12004 


2014D have shown that 


this second kind of pattern can be explained by the su¬ 
perposition of two components: an H-component pro¬ 
ducing the heavy r-process elements (and maybe also 
lighter ones) and an L-component contributing only be¬ 
low Z = 50. In neutron star mergers, the dynamic and 
viscous ejecta can account for the H-component while 
the wind ejecta would be the L-component. Matter with 
a perfect mixing of the three ejecta will lead to a pat¬ 
tern like in CS 22829-052. This is shown in the upper 
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panel of Fig. where the combination of wind (L- 
component) ana dynamic (H-component) ejecta repro¬ 
duce the observed abundances. The differences around 
the thi rd peak are due to n uclear physics input as shown 
~ ' (2015). If the mixing is not perfect 


Eichler et al. 


m e.g. 

and the wind ejecta combines only with a small amount 
of dynamic and/or viscous ejecta, then one can explain 


the tre nd observed in HD 122563 (see also |Just et ah 
(|2015a|)). This is shown in the bottom panel of l^lg. |l2 
where the dynamic ejecta has been reduced by a factor m 
50. In this case, one would expect to have more variabil¬ 
ity in the observations depending on the amount of the 
H-component that is mixed. Indeed observat ions show 


variability for Z > 50 in these kind of stars (jRoederer 
et al.||201Q Hansen et al.||2014 ) 



z 



Fig. 12.— Top: combined elemental abundances of neutrino- 
driven wind and dynamic eject a in comparison to t he yields of the 
metal-poor star CS22892-052 (|Sneden et al.||2003|). Abundances 
are weighted with the ejecta mass to consistently add the two con¬ 
tributions. Bottom: c omparison of the combined yields with the 
ones of the HD stars (|Wallerstein et al.||1963| |Honda et al.||2004| 

a . Note that the dynamic ejecta are scaled down by a lactor 


On the time scale of star formation, it is likely that 
all the ejecta have efficiently mixed, even if the wind 
expands perpendicular to the disk and the other two 
are more homogeneously distributed or stay in the disk 
plane. Do we expect stars to be polluted by a neu¬ 
tron star merger preferentially only receiving axial ejecta 


rather than disk ejecta? Or do Honda-like abundances 
(L-component) come from different sources (e.g. core¬ 
collapse supernovae with slightly neutron-rich condi¬ 
tions)? More investigation of the long-term mixing is 
necessary to understand the contribution of neutron star 
mergers to the different patterns found in the oldest ob¬ 
served stars. However, the dynamic and wind ejecta will 
not mix before they are substantially decelerated by the 
dilute ambient medium. For example, assuming a par¬ 
ticle density of ^ I cm“^ for the ambient medium, the 
estimated deceleration time scale is of the order of several 
years. The dynamic ejecta expand earlier and with high 
initial velocities of ^ O.Ic. The neutrino-driven wind, 
in contrast, gets unbound after neutrinos have deposited 
enough energy. Therefore, the wind matter is ejected 
later and eventually reaches asymptotic velocities of up 
to ^ 0.08c. The different ejection times and velocities 
make it possible to observe a kilonova signal from differ¬ 
ent parts of the ejecta with different composition, as we 
discuss in the next section. 

3.3. The electromagnetic signal: a semi-analytical model 

The nucleosynthesis network allows us to compute ra¬ 
dioactive heating rates for the wind outflow. Fig. [T^ 
shows the heating rates for different bins, normalized to 


en(0 = er 

2 ; g ^s ^ with time tc 

] in units of 

days (jMetzger et al. 

2010 

Korobkin et al.| 

2012 

). It is 
aeating 

interesting to point out that all the normaJ 

iized 


rates show considerable excess at t ^ 4 h. This excess 
is caused by decay of radioactive nuclei in the vicinity 
of the first r-process abundance peak {A ^ 80), as those 
possess the highest mass fractions. The properties of the 
isotopes which contribute most to the heating rates in 
the neutrino-driven wind are listed in Tab. [TJ On aver¬ 
age, about 40% of the decay energy is carried away only 
by neutrinos. An additional fraction of energy is radiated 
away by escaping photons, while both the electrons and 
the remaining part of the photons thermalize. Therefore 
our hypothesis of an effective beta-decay thermalization 
fraction of 50% is reasonable. 


TABLE 1 

Properties of the dominant /3-decay nuclei at t ~ 1 day; 

BASED ON THE DATA FROM NuDAT 2.6 DATABASE 

(http://www.nndc.bnl.gov/nudat2/) 


Isotope 

*1/2 

(h) 

Qia) 

(MeV) 




T^avg (e) 

(MeV) 

88 Kr 

2.83 

2.92 

0.12 

0.21 

0.67 

1.34 

88 Rb 

0.30 

5.31 

0.39 

0.49 

0.13 

1.59 

8'^Kr 

1.27 

3.89 

0.34 

0.46 

0.20 

0.95 

83Br 

2.37 

0.98 

0.33 

0.66 

0.007 

0.46 

81 Sr 

0.37 

3.93 

0.28 

0.37 

0.35 

0.42 

i'8Ge 

1.47 

0.96 

0.24 

0.47 

0.29 

0.28 

1'8As 

1.51 

4.21 

0.30 

0.39 

0.31 

0.94 

1’1’Ge 

11.2 

2.70 

0.23 

0.36 

0.41 

0.47 


(a) Total energy released in the decay; (b),(c),(d) Fraction of 
the decay energy released in electrons, neutrinos, and 7—rays; 
(e) Average photon energy produced in the decay. 


Her e we e xtend the semi-ana l ytic in odel of |Grossman| 
et al. (2014) and Perego et al. (2014) to conipute light 


curves. We assunie that the wind shuts off immediately 
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Cylindrical Radius R [10^ km] 


Fig. 13.— Radioactive heating rate, calculated with the nu¬ 
cleosynthesis network for a representative tracer of the neutrino- 
driven wind in each angular bin, normal ized to the fitted expression 


eo{t) = e rg g ^s ^ (see, e.g., Metzger et al. ( 2010| >; 


|robkin et al 
by a tact or 


Ko- 


T|(|2012|)). In all cases the heating rate exhibits an excess 
or ~ 2.5 at t ~ 4 h. 


Fig. 14. — Distribution of summed up mass fractions for nuclei 
with mass number A > 130 in the neutrino-driven wind. The 
upper half space is shown in cylindrical coordinates for a time of 
50 s after the merger. Tracers with vanishing total mass fractions 
Xheavy < 10“^° are set to 10“^°. 


after the collapse of the MNS and we explore collapse 
times of 90 ms, 140 ms and 190 ms after merger. The 
masses as well as average properties of the ejecta in the 
four bins for each of the cases are given in Tab. The 
higher value of /^ = 10 cm^/g for the opacity in the last 
angular bin is justified by inspecting the distribution of 
mass fractions for heavy elements Xheavy (foi* mass num¬ 
bers A > 130) in the neutrino-driven wind. This dis¬ 
tribution is portrayed for extrapolated positions of the 
unbound tracers in Fig. [^at t = 50 s after the merger, 
saw el( 


Mass fractions of heavy elements are marginal in the 
bins 1 — 3, but heap up in the area of bin 4. Reflect¬ 
ing that sizeable fractions of heavy elements are present, 
we adopt an opacity similar to the one of lanthanides or 
actinides ( Kasen et al.||2013 ) for the last bin. 


TABLE 2 

Parameters used for computing luminosities of individual 

BINS: SPANNED SOLID ANGLE, MASS OF THE BIN, AVERAGE 
ASYMPTOTIC VELOCITY OF THE EJECTA IN THE BIN, AND THE 
ADOPTED VALUE OF THE GRAY OPACITY. 


^sim 

bin 

AQ 

Amej/10 ^Mq 

'^e.i/c 

R [cm^/g] 

90 ms 

1 

0.21409 

0.013 

0.05055 

1 


2 

0.62769 

0.30 

0.07974 

1 


3 

0.99851 

0.77 

0.07287 

1 


4 

1.30129 

0.34 

0.06808 

10 

140 ms 

1 

0.21409 

0.05 

0.04655 

1 


2 

0.62769 

0.75 

0.07483 

1 


3 

0.99851 

1.99 

0.07626 

1 


4 

1.30129 

2.43 

0.06694 

10 

190 ms 

1 

0.21409 

0.07 

0.04694 

1 


2 

0.62769 

1.18 

0.07325 

1 


3 

0.99851 

3.39 

0.07500 

1 


4 

1.30129 

4.80 

0.06466 

10 


Figure [15] illustrates our semianalytic model for light 
cu rve calculations, combi ning the results of this work and 
of Rosswog et al. (2014). The neutrino-driven wind is 
schematically shown in blue, while the gray-shaded areas 
tag the density of the dynamic ejecta in steps of 0.5 dex. 


We subdivide the wind outflow into the same four bins 
as in Section 2.1 above. Each bin is then approximated 
by a conical slice of a spherically-symmetric outflow with 
radial density distribution averaged over the bin and ex¬ 
panding wit h averaged velocity (a s imilar approach is 


employed inlM 
dio flares). AT 


jMargalit & Piran (2015) in a context of ra- 
^11 radiation is assumed to escape from the 
photosphere for each bin. Vector ft is the unit normal to 
the photosphere, and an observer is pointed to by a unit 
vector q. The sketch also illustrates possible obscuration 
effects from the dynamic, very opaque ejecta, which are 
concentrated in a puffed up toroid around the equatorial 
plane. It is apparent that when viewed from the equa¬ 
torial plane, the wind emission is completely obscured, 
while from the pole it is possible to observe at least the 
first three bins. 



Fig. 15.— Morphology of the neutr ino-driven wind fblue. this 
work) and dynamic ejecta (gray, cf. [Rosswog et al. ijM) in 
asymptotic velocity space, represented by a polar cut througH the 
outflows. Shaded areas represent the density of dynamic ejecta, 
spaced by 0.5 dex from the maximum density at about vji ~ 0.115c. 
Each angular bin is approximated by a conical slab with unit vector 
fi normal to the photosphere. An observer is pointed to by a unit 
vector q. 
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For a bin k with mass Arrik which spans solid an¬ 
gle 2A11/C (a factor of 2 takes into account upper and 
lower lobes), we take a spherically-symmetric outflow 
with mass An Arrik/2Aftk and compute i ts isotropic bolo- 
metric luminosity Lk iso as described in Grossman et ah 
dMil ) (see their Sect. 4.1) 


I'he luminosity is generated by radioactive heating in 
the bulk of the outflow and the resulting photons escape 
from a photosphere. Spanning the same solid angle 2A^k 
for each bin, the luminosity for bin /c is a proportional 
fraction of isotropic luminosity of the spherical model: 

j _ AQ^k j 

We then compute the total luminosity of the wind out¬ 
flow by summing up individual bin contributions, ignor¬ 
ing possible radiative flux between the bins. Individual 
bin contributions for the combined bolometric luminos¬ 
ity of the neutrino-driven wind and dynamic ejecta for 
three cases of the MNS collapse time delay are displayed 
in Fig. For the dyn amic ejecta, we used an aver¬ 
age case(model A from Grossman et al. ( 2Q14[ )) with 
mass 1.3 • 10“^ Mq. Notice that the luminosity of bin 1 
is dimmer than for the bins 2 and 3, because it spans 
a smaller solid angle. The contribution from bin 4 is 
not only smaller, but also peaks much later (3 — 4 days) 
due to the hi gh opacity caused by the very heavy nu¬ 
clear content (Kasen et al.||2013), such that the medium 
becomes delayed transparent when the temperature has 
decreased. 



0.01 0.1 1 0.01 0.1 1 10 


time [d] 


Fig. 16.— Individual contributions of the combined ejecta to 
the luminosity in the angular bins for three different MNS collapse 
times: 90 ms (solid), 140 ms (dashed) and 190 ms (dotted line). 

Combined bolometric luminosities of the neutrino- 
driven wind and dynamic ejecta are shown in Fig. EZl 
A first peak of Apeak ^ 4 • 10^^ erg/s is reached at about 
4 h after the merger. This luminosity stays roughly con¬ 
stant for several days while rapidly shifting into the in¬ 
frared band, as shown below. For higher MNS collapse 
times, the light curve exhibits a double peak structure, 
but overall the bolometric light curves on the plots ex¬ 
hibit no appreciable difference between different MNS 
collapse times. 

Our simple model also allows to calculate approximate 
light curves in different bands for different orientations 



Fig. 17.— Bolometric luminosities for three cases of MNS collapse 
times: 90 ms (solid), 140 ms (dashed) and 190 ms (dotted line), 
computed by summing up contributions from individual bins and 
adding the contribution from dynamic ejecta. Contribution from 
the latter is also plotted separately (dash-dotted line). 


of the system with respect to an observer. At any given 
moment t we can compute an effective temperature T/^^eff 
for a bin k using the formula: 

Lk,iso{t) = 47rr| ph ■ (crsBT^.eff) (3) 

where rk,ph{t) is the photospheric radius for the spher¬ 
ical model of bin k, and ctsb is the Stefan-Boltzmann 
constant. We make an assumption that all the flux is 
emitted from the photosphere with blackbody spectrum 
Bu{Tk^efi) and isotropic intensity (following Lambert’s 
law). Then, the spectral flux in the direction q is given 
by: 

Mq, i) = E MTkMt)) [[ (9 ■ dn), (4) 

k ^ (^ fc -^>0 


where fik is unit normal to the photosphere of bin k, and 
integration for each bin spans the part of the surface fac¬ 
ing the observer (stated by the condition {fik * ^ > 0). 
This is also illustrated in Fig. Essentially, the inte¬ 
grals in this formula are simply time-independent geo¬ 
metric projections of the photosphere for each bin onto 
the viewing plane. They can be computed beforehand 
and used as weighting factors for calculating the 

combined light curve: 

= 'Y^Pk{q)Bu{Tk,es{t))- (5) 

k 


After integrating over certain frequency ranges, 

we obtain broadband light curves, shown in Fig.[^ The 
left-hand panel shows synthetic light curves for wind out¬ 
flow only, while the right-hand panel displays the com¬ 
bined contribution from both wind and dynamic ejecta. 

Orientation effects are displayed as a range of values 
for each band, with the maximum magnitude achieved 
when the system is viewed ”face-on”, i.e . from the pole. 
As pointed out by Kasen et al. (2015), even a small 
amount of dynamic ejecta can completely obscure the 
optical emission from the winds. Therefore, when adding 
up contributions from different bins and dynamic ejecta, 
we emulate this obscuration by excluding certain bins, 
depending on the inclination angle 0 with respect to the 
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observer. Our choice of which bins to exclude is moti¬ 
vated by the geometry of the ejecta (depicted schemat¬ 
ically in Fig. 15). Specifically, we assume that bin 4 is 
completely obscured for all observing angles; bins 2 and 
3 are obscured when 0 > 60°, and for the edge-on view 
with 0 = 90° the wind outflow emission is obscured com¬ 
pletely. For this reason, depending on system orienta¬ 
tion, the luminosity in different bands can vary by up to 
an order of magnitude, from bright blue when observed 
from the pole, to dim infrared when observed from the 
side. 



0 5 10 0 5 10 

time [d] 


Fig. 18.— Broadband light curves of the wind outflow (left panel) 
and wind-hdynamic ejecta (right panel), showing the combined blue 
(U-hV bands), red (V-hR) and infrared (J-hH-hK) contributions. 
The top, middle and bottom rows show the three cases of MNS 
collapse times: 90 ms, 140 ms and 190 ms, respectively. The range 
for each light curve spans possible inclination angles of the system 
axis with respect to the observer: 0° (solid) - view from the top, 
30° (dashed), 60° (short dashes) and 90° (dotted) - view ’’edge- 
on” . 


4. DISCUSSION AND CONCLUSIONS 

We have presented a comprehensive study of the 
neutrino-driven wind nucleosynthesis in the aftermath of 
a double neutron star merger. We have focused on the 
initial phase where the remnant consists of a MNS sur¬ 
rounded by a massive accretion disk. The current nucle¬ 
osynthesis study is based on the first three-d imensional 
simulation of such a wind (Perego et al. IMit , sampling 
the wind outflow by ^ 17UUU tracer particles. For ev¬ 
ery tracer, nucleosynthesis calculations are performed in 
a post-proces sing step. This wor k advances the prelimi¬ 
nary study in|Perego et al.|(|2014|) that only used 10 trac¬ 
ers and thus led to a slight overestimation of the amount 
of heavy elements (^4 > 130) produced, whereas the total 
mass ejected in the wind was underestimated. Our main 
results are summarized and discussed alongside their un¬ 
certainties in the following. 


The abundance distribution in the neutrino-driven 
wind is characterized by the production of lighter heavy 
elements (^4 < 130). We have found a time and angle 
dependency of the composition correlated with the elec¬ 
tron fraction evolution. Elements heavier than A > 130 
are produced only at early times, in an almost negligible 
amount and in the wind region that is closest to the disk 
(between 45° and 60° from the pole, bin 4 in Fig. 0. In 
contrast, matter ejected more perpendicular to the disk 
(bin 1: 0° — 15°) produces only nuclei up to the first r- 
process peak (^4 < 100). The nucleosynthesis path of the 
neutrino-driven wind after a neutron star merger does 
not reach the heaviest nuclei, in some cases it even stops 
before N = 82. Relevant nuclear physics input (e.g., nu¬ 
clear masses and beta-decays) along the path is known 
or will be measured in current and future generation fa¬ 
cilities in the coming years to provide more reliable data. 
Our calculations ar e based on the FRDM mass model 
(Moller et al. 1995|) and consistently calculated nuclear 


reaction rates ([Rauscher fc Thielemann 2000). A de- 


can g 

:ive rise to new insl 

ghts (de Jesus Mendoza-Temis 

et al.| 

2014 Eichler et al.( 

|2015 

Goriely et al.l|2015). Es- 


and neutron captures may have an impact on the abun¬ 
dances before the second r-process peak where we find a 
significant trough. 

The electron fraction is here the key parameter to un¬ 
derstand the characteristics of abundances. It covers a 
broad range from 0.1 to 0.4 with average values of ^ 
0.33. As noted already in earlie r studies (see in particular 
Fig. 8 in Korobkin et al!](|2012 ) where the Te-dependence 


was e xplored; |Freiburghaus et al.||1999a| |Kasen et al.| 


2015), Yq « 0.25 is the threshold beyond which no 


more heavy r-process elements are formed. Since abun¬ 
dances depend strongly on 1^, an accurate treatment of 
weak reactions is crucial and can be achieved by improv¬ 
ing the neutrino transport and including more neutrino- 
matter reactions. We have estimated a ^ 20% error on 
the electron fraction due to the approximate transport 
scheme. This is suggested by the compa rison of neu¬ 
trino l umin osities and mean ene rgies from [Perego et ah] 


(2014) and Dessart et al 
best treatments tor these 


(2009), which are amon g the 


_ kinds of systems (se e also Just 

et al. ( 2015a|b ); Fernandez fc Metzgerj ( 201^ ). For the 


neutrino-matter interactions, a small variation of can 
arise from a detailed inclusion of the (anti)neutrino ab- 
sorption including weak magnetism and medium effects 


(Roberts et al.||2012 Martmez-Pinedo et al. 

T'he amount of ejecta also depends on angle and time. 
As the wind is preferentially driven out at lower latitudes, 
bin 3 and 4 contain roughly 6 times the matter of bins 1 
and 2. About 110 ms after the merger a steady-state is 
reached in terms of mass loss. The total unbound mass 
depends crucially on the time when the MNS collapses 
and if it survives for ^ 200 ms as much as 9 • 10 
can become unbound (i.e. ^ 5% of the initial mass of 
the accretion disk). This amount of matter rivals the 
dynamic ejecta that are typically found to be of order 
of O.OIMq (1.3 • 10~^ M^ for the specific case considered 
here). Recent studies ([Fbrnandez fc Metzger|2013||K asen 
et al.||2015l|Just et al.||2015a[ Sekiguchi et a[.||2015[) indi- 


cate that the winds switches oft rapidly once a black hole 
forms. Approximately half of the neutrino luminosity is 
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provided by the MNS. Therefore, when the black hole 
forms this contribution vanishes and the disk structure 
may significantly change. It is possible that general rel¬ 
ativistic effects also have a non-negligible impact on the 
disk winds from a MNS and accretion disk sy stem, as is 


the case for neutron star - black hole systems ([Caballero 
eFaI][M^ . 


I'he angle dependency of the composition may have 
consequences for the mixing with other ejecta. Both dy¬ 
namic and viscous ejecta contribute also to the produc¬ 
tion of heavy r-process elements beyond the second peak 
(^4 > 130). If wind material perfectly mixes with these 
ejection channels, its contribution may lead to variations 
for the lighter heavy elements (i.e., A < 130). This vari¬ 
ation is indeed observed in stars with high enrichment of 
heavy r-process elements. However, if the mixing is not 
perfect one could speculate that neutron star mergers 
with long-lived MNS may produce different abundance 
patterns including one with low enrichment of heavy r- 
process el ements. This kind of patt ern (“Honda-like pat¬ 


tern”, see Honda et al. (2004 2007)) is observed in a few 
very old stars and its origin is still unknown. However, 
detailed mixing models following the late evolution of the 
three nucleosynthesis-relevant ejecta are required before 
concluding about the contribution of neutron star merg¬ 
ers to observed abundance patterns. 

The initial separation of wind and dynamic ejecta 
allows to describe the individual contributions to the 
light curve of the electromagnetic transient. We use a 
semi-analytic radiation transport model to calculate the 
luminosities of the wind and dynamic ejecta. The light 
curve of the wind ejecta peaks after ^ 4 hours in the 
blue compared to the one of the dynamic ejecta that 
peaks few days after the merger in the infrared. The 
combined light curve will significantly depend on the 
observing an gle, as previously pointed out by [Kasenj 


et al. 
of laiit' 


([20T^. 

thaniae 


Due to the hig h line expansion op acities 
es and actinides (Kasen et al. 2013), even 


small amounts of dynamic ejecta along the line of 
sight can completely obscure the blue transient from 
the wind. Since the morphology of dynamic ejecta is 
such that it leaves polar regions evacuated, it is far 
more likely to detect the blue component when the 
system is observed from the pole. In Fig. 18, we have 
presented the light curves of obscured and unobscured 
winds as a function of angle and time of black hole 
formation. The luminosity in the blue band exhibits a 
much higher dependency on the viewing angle compared 
to the luminosit y in the infrared, i n agreement with 
previous studies (Kasen et al. 2015). This significant 
anisotropy needs to be taken into account when assessing 
detectability of optical counterparts to neutron star 
mergers. Another interesting feature is the noticeable 
dependence of the position of optical peaks on the 
time of black hole formation (apparent in Fig. 18). A 
longer collapse time leads to a later and more prolonged 
peak in both U+B and V+R bands, while this trend is 
washed out in the infrared. Although at this point it 
is difficult to draw any quantitative conclusions due to 
limitations of our model, our results suggest a possibility 


^ Despite these differences, the overall amount of dynamic ejecta 
is in good agreement between different studies. There is a depen¬ 
dence on the mass ratio and the equation of state, with softer EOSs 
ejecting more mass in GR calculations, but ejecta masses up to a 
few percent of the solar mass can be reached, see, e.g. Tab. 1 in 


to use the position and duration of optical transients as 
proxies for the time of blac k hole formation (see also 
Metzger & Fernandez (2014|)). More work is needed to 
take into account additional channels of late-time ejecta 
pro duction, and to im prove radiation treatment (see. 


e.g. 


Kasen et al. 


(2015)). 


Our results clearly indicate the importance of neutrino- 
driven winds for the nucleosynthesis of neutron star 
mergers. It needs to be stated, however, that these re¬ 
sults have been obtained based on Newtonian hydrody¬ 
namics simulations and that general-relativistic effects 
could potentially modify some of ou r conclusions. For 
example, recent studies in full G R (Hotokezaka et al. 
2011, 2013 Sekiguchi e t al. 2015|) ana m the conformal 
flatness approximation ( [Banswein et al.|2013 ) agree that 


shocks are more important lor the dynamic ejecta in the 
GR case, especially when a soft equation of state is used. 
This is mainly because the stars are - for a given EOS - 
more compact due to the stronger relativistic gravity and 
hence lead to larger velocities prior to contact. A softer 
EOS leads compared to a stiff one to more compact stars 
and also to smaller sound speeds which, in turn, make 
the occurrence of shocks more likely. Therefore, the dy¬ 
namic ejecta are dominated by the hotter ’’interaction 
component” while Newtonian ejecta are dominated by 
the co lder, unshocked ’’tidal component” (Korobkin et al. 
2012) ^ Therefore, a GR calculation will plausibly cause 
a higher temperature environment which is expected to 
lead to larger electron fractions in the ejecta. Due to 
shock-heating it is also likelier that unbound material is 
ejected out-of-plane. Thus a neutron star merger could 
actually, consistent with our findings he re, produce a 
very broad range of r-p rocess abundances (Wanajo et al 
20141 [Just et al.||2015a|). 

Dnderstanding the interplay of all ejecta from neutron 
star mergers is key to decode the diverse aspects of this 
astrophysical site. Eurther detailed simulations and com¬ 
plete nucleosynthesis studies will help to pin down the 
role of neutron star mergers for the origin of the heavy 
elements in the universe as well as to predict reliable 
electromagnetic counterparts for the detection of GWs. 
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